Lasing and transport in a coupled quantum dot-resonator system 
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We investigate a double quantum dot coupled to a transmission line resonator. By driving a 
current through the double dot, a population inversion between the dot levels can be created, and 
a lasing state of the radiation field is generated within a sharp resonance window. The transport 
current correlates with the lasing state. The sharp resonance condition allows for resolving small 
differences in the dot properties. Dissipative processes in the quantum dots and their effect on the 
lasing and transport behavior are investigated. 

PACS numbers: 42.50.Pq, 73. 21. La, 03.67.Lx 
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I. INTRODUCTION 

The interaction between light and matter is one of the 
fundamental topics in physics. Modern techniques al- 
low for fabricating artificial atoms with highly control- 
lable parameters, which provide an ideal testing ground 
for quantum optics. In the emerging field of circuit 
quantum electrodynamics (QED) [lH5[, where super- 
conducting qubits, playing the role of artificial atoms, 
are strongly coupled to a superconducting transmission 
line resonator, many quantum optics effects have been 
demonstrated with unprecedented accuracy, and some 
novel phenomena have been observed. An example is 
single-qubit lasing [6 8] , where quantum noise influences 
the linewidth of the emission spectrum in a characteris- 
tic way [§-[l2j]. Recently proposals to replace in circuit 
QED setups the superconducting qubits by semiconduc- 
tor quantum dots were put forward (13l - Fl7l j . Quantum 
dots have the advantages that many of the parameters, 
e.g., the level structure, coupling strength, and tunneling 
rates, can be varied in situ, and details of the transport 
properties can be resolved via adjacent quantum point 
contacts [l8l - f20l |. Fast experimental progress has been 
made towards coupling quantum dots to a superconduct- 
ing resonator [2l| - r23j . 

In this paper we consider such a circuit QED setup 
with a semiconductor double quantum dot coupled to a 
high-Q transmission line. We find that by driving a cur- 
rent through a suitably biased double dot, a population 
inversion can be created between the dot levels, which 



then may lead to a lasing state in the resonator. This 
lasing state in the resonator correlates with the trans- 
port properties through the double dot system. This 
property not only allows probing the lasing state via a 
current measurement, which may be easier to perform 
in an experiment, but also, because of the narrow res- 
onance window, opens perspectives for applications for 
high resolution measurements. 

The paper is organized as follows. In Sec.|TI]we present 
the model for the coupled quantum dot-resonator system 
and the pumping mechanism, provide estimates for the 
parameters, and we introduce the theoretical tools for 
the analysis of the system. In Sec. IIIII we investigate 
the stationary properties of the radiation field inside the 
resonator, while its emission spectrum and linewidth are 
discussed in Sec. [IVJ In Sec. [V] we analyze the correla- 
tion between the lasing and transport properties of this 
coupled system. We conclude with a brief summary. 



II. MODEL AND METHOD 

We consider a system shown schematically in Fig. [TJ 
The double quantum dot is biased such that three charge 
states play a role, a reference neutral state and two states 
differing by a single extra electron in the left or right dot. 
The three states are denoted by |0, 0), |1,0), and |0,1), 
respectively. The two singly charged states, with energy 
difference e (denoted as dot detuning), are coupled by 
coherent interdot tunneling with strength t. The double 
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transmission line oscillator 




FIG. 1: (Color online) Coupled quantum dot-resonator cir- 
cuit. The double dot is placed at a maximum of the electric 
field of the transmission line in order to maximize the dipole 
interaction with the resonator. For gate-defined quantum dots 
most parameters can be tuned by the applied voltages. 

dot system is therefore described by a two-level Hamil- 
tonian, 



(1) 



with Pauli matrices t z = |1, 0) (1, 0| — |0, 1)(0, 1| and t x — 
|1,0)(0,1| + |0,1)(1,0|. 

We model the superconducting transmission line as a 
harmonic oscillator with frequency tj T> and a and cv de- 
noting the annihilation and creation operators for the ra- 
diation field. In typical circuit QED experiments, the res- 
onator frequency is of the order of several GHz, and the 
vacuum-fluctuations-induced voltage drop V T between its 
central strip and ground plane is of the order of pN pj . 
For definiteness we choose lu t /2tt — 4 GHz throughout 
this paper. For the geometry shown in Fig. [TJ the su- 
perconducting resonator is capacitively coupled to the 
charge states of the double dot which have different 
dipole moments. It leads to the interaction 



H c = hg Q (a f + o)t z 



(2) 



When the interdot tunneling is weak and the single- 
particle wavefunctions are strongly localized in either 
dot, the coupling strength can be estimated as go ~ 
eEd/(2h) ~ 50 MHz, where d ~ 0.3 pm is the distance 
between the centers of the two dots and E ~ 0.2 V/m the 
electric field at the antinode of the resonator mode. This 
estimate agrees with the experimental observations [24] ]. 

In the eigenbasis of the two- level Hamiltonian (p}, 
within the rotating wave approximation, the Hamilto- 
nian for the coupled system reduces to 



a z + Hu> T a' a + ftg(a^er_ + aa + ). 



(3) 



Here a z = |e)(e| — \g)(g\ (similar for a±) is defined in this 
subspace spanned by eigenstates 



|e) = cos (0/2)|l,O)+sin(0/2)|O,l), 
\g) = -sin(0/2)|l,O) + cos(0/2)|O,l) 



H EL = -eV 
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FIG. 2: (Color online) Tunneling sequence in the double dot 
system. The chemical potentials /xl(1,0) and pn{0, 1) of the 
two dots are assumed to be arranged as indicated. 



with angle 9 = arctan(£/e) characterizing the mixture 
of the charge states. The energy splitting between the 
eigenstates, Hloq = ^/e 2 + t 2 , can be tuned via gate volt- 
ages, allowing control of the detuning relative to the 
resonator frequency A = u)q — w r . The effective dot- 
resonator coupling strength, g = g Q sin0, reaches its max- 
imum at e = 0. To achieve a lasing state, the detuning 
A should be small and the coupling g strong enough to 
overcome the dissipative processes. This requires both 
the energy difference e and the interdot tunneling t to be 
of the order of several /xeV. 

In order to create a population inversion, which is usu- 
ally required for lasing behavior, we consider a pumping 
scheme involving the third state |0, 0), similar to an op- 
tical laser. (Note that lasing without inversion is also 
possible in circuit QED setups where the dissipative en- 
vironment enhances the photon emission as compared to 
absorption 25]. Population inversion induced by other 
pumping mechanisms, e.g., driving the two- level system 
coherently @, H3) were also discussed.) As indicated in 
Fig. [2l we assume that a bias voltage is applied across 
the double dot such that only the chemical potentials of 
the states 1 1 , 0) and |0, 1) lie within the bias window. At 
low temperature compared to the charging energy, the 
only possibility for an electron to tunnel into the dot 
system is from the left lead to the left dot, leading to the 
transition from state |0, 0) to |1,0) with tunneling rate 
IY,. Similarly, an electron can tunnel out into the right 
lead, creating a transition from |0, 1) to |0, 0) with tun- 
neling rate Tr. The pumping leads to a non-equilibrium 
state where the population in the state |1, 0) with higher 
energy is enhanced compared to that in |0, 1). The pop- 
ulation inversion persists when we go to the eigenbasis, 
where it is given by tq = 4cos6*/[3 + 008(26*)] fl71 ]. 

To analyze the dynamics of the coupled dot-resonator 
system, we employ a master equation for the reduced 
density matrix p in the Born-Markovian approximation, 



i 

Aot P ■ 



(5) 



(4) 



We assume that the dissipative environment is charac- 
terized by smooth spectral functions [13, Hll and can be 
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described by Lindblad operators of the form 
C lP =^(2L lP L\-L\L lP -pL\L, 



(6) 



For the oscillator we take the standard decay terms 
L r = a with rate T r = n. Throughout this paper we con- 
sider low temperatures, T = 0, with vanishing thermal 
photon number and excitation rates. For the two-level 
system we account for the relaxation by L± = <r_ with 
rate T± and for pure dephasing by L v — a z with rate F* . 
For the discussion of lasing behavior, it is worthwhile to 
introduce a total decoherence rate for the dot system, 
namely, T v = T i /2 + T sin 4 (<9/2) + cos 4 (<9/2) /2 + r* . 

The last two terms account for the incoherent tunneling 
between the electrodes and the left and right dots. In the 
basis of the pure charge states, the Lindblad operators 
describing these processes are given by £l = |1,0)(0, 0| 
and £r, = |0, 0)(0, 1|. For simplicity, we set Tl = Tr = T 
throughout the paper. 



III. STATIONARY PROPERTIES 

We start with analyzing the stationary properties of 
the radiation field inside the resonator, namely the av- 
erage photon number (n) and the Fano factor F = 



")/(n), which serve as good indicators of the 



state of the radiation field. We plot both as functions of 
the detuning in Fig. [3] To begin with we neglect the re- 
laxation and decoherence of the dot levels, = T* = 0. 
When the detuning is large (e.g., the point/curve marked 
by 'a' in Fig. [3]), the quantum dot effectively does not in- 
teract with the resonator, and the photon number is low 
with a thermal distribution. For weaker detuning (e.g., 
point 'b') the photon number increases sharply and sat- 
urates at resonance, A = 0. Approximately, the photon 
number can be estimated as [29(, 



F cos cos(26») + 7 



3 K 



96 g 2 



(4 A 2 + r 2 ). 



(7) 



In the typical parameter regime of experiments the peak 
value of the photon number is mainly given by the first 
term, i.e., the ratio between the incoherent pumping and 
the decay rate of the resonator. For large photon number 
good agreement is reached between the analytical expres- 
sion (dotted line) and numerical results (solid line). 

Statistical properties of the radiation field can be ob- 
tained by investigating the Fano factor. When the res- 
onator is in a lasing state (e.g., point 'c'), the radia- 
tion field is in the coherent state and the photon num- 
ber distribution is Poissonian. In this case the Fano fac- 
tor equals to 1. In the strong-coupling regime, the Fano 
factor can even become smaller than I (e.g., point 'd'), 
which indicates a sub-Poissonian distribution of the radi- 
ation field. In this non-classical regime, the photon num- 
ber distribution is squeezed compared to the Poissonian 
distribution. 
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FIG. 3: (Color online) Upper panel: Average photon num- 
ber (n) and Fano factor F as functions of dot detuning e for 
vanishing decoherence of the dot levels (r^ = TZ, = 0). The 
analytical approximation (n) a (dotted line) in Eq. (J7J is com- 
pared to the numerical result (solid line). The decay rate of 
the resonator is chosen to be k/2tt = 40 kHz. Throughout this 
paper we choose the incoherent tunneling rate V = 4MHz, the 
coherent interdot tunneling strength t — 5 /xeV, and the bare 
coupling strength go = 4MHz. Lower panel: Photon number 
distribution at points a, b, c, and d marked in the upper panel. 
Dotted line indicates a Poissonian distribution for comparison 
with the sub-Poissonian distribution at point d. 



We now turn to the effects of decoherence of the dot 
levels on the lasing state, which are illustrated in Figs. [4] 
The relaxation process reduces the pumping efficiency 
and deteriorates the lasing state. Hence the photon num- 
ber decreases. At the same time the Fano factor around 
resonance becomes larger than I, resulting from an in- 
crease of amplitude fluctuations. Pure dephasing, on one 
hand, diminishes the efficiency of energy exchange be- 
tween the dot levels and the resonator. On the other 
hand, it effectively broadens the window in which the 
dot levels can interact with the resonator. As a result, 
the peak of the photon number is broadened, while the 
peak value decreases. 
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FIG. 4: (Color online) Average photon number (n) and Fano 
factor F as functions of detuning. Solid lines represent the 
results without decoherence, blue dotted lines are those for 
relaxation rate Fi/2n = 1.2 MHz and vanishing pure de- 
phasing, while red dashed lines are for pure dephasing rate 
F*/27r = 12 MHz but vanishing relaxation. Other parame- 
ters are chosen to be the same in Fig. [3] 

IV. EMISSION SPECTRUM AND LINEWIDTH 

We now turn to the spectral properties of the radiation 
field. The emission spectrum of the radiation field, 

/oo 
e-*"V(t)a(0)), (8) 
-oo 

which can be measured directly in experiments, provides 
further insights into the state of the radiation field. To 
calculate the phase correlator, we make use of the quan- 
tum regression theorem [27j , 

(ot(t)o(0)) =Tr[a?e c ^ t ap st ], (9) 

with p st being the stationary solution of the master equa- 
tion (g§. 

The spectral properties of the radiation field in the 
quantum dot-resonator system are similar to those of a 
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FIG. 5: (Color online) Emission spectrum S e (tj) (in arbitrary 
units) for different detunings. 
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FIG. 6: (Color online) Linewidth of the emission spectrum 
as a function of detuning. The green solid line is the result for 
vanishing relaxation and decoherence Fj. = FJ, = 0, the blue 
dotted line for relaxation rate F^/27r = 1.2 MHz but vanishing 
pure dephasing, and the red dashed line for pure dephasing 
rate F*/27r = 4 MHz but Fj. = 0. Here the decay rate of the 
resonator is chosen to be k = 80 kHz. 

single-qubit maser discussed in Refs. [Ml. In Fig.[S] 
we plot the emission spectrum for different values of the 
detuning. At resonance, A = 0, the emission spectrum 
shows a sharp peak at the resonator frequency io T . For 
nonvanishing detuning the peak of the emission spectrum 
is shifted in frequency by 8u). For weak detuning, this 
shift grows linearly with the detuning, Suj ~ An/2T V , 
with T v =rj,/2 + r/4 + T* characterizing the total dis- 
sipation rate of the quantum dot system. At the same 
time, the peak value of the emission spectrum decreases 
with growing detuning, related to a lowering of the aver- 
age photon number. 

The linewidth of the emission spectrum is plotted as 
a function of detuning in Fig. [6] In a typical situation 
where the coupling g is not too strong compared to the 
decoherence of the coupled system (k, the linewidth 
grows monotonously with increasing detuning, as indi- 
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FIG. 7: (Color online) Current as a function of dot detuning 
e. The parameters are chosen to be the same as in Fig. O The 
dot levels become resonant with the resonator at 15.9 fieV. 



cated, e.g., by the blue dotted line. At resonance, deep 
in the lasing regime, the linewidth can be estimated as 



k l oc cry (r ¥ 



i) m 



It becomes very small for large 



photon numbers. Here, due to the restrictions of our 
numerical calculation, we illustrate to a situation with 
low photon number around 15 and not very small val- 
ues of «l- For stronger detuning the system leaves the 
lasing regime, and the linewidth is approximately given 
by kl — — g 2 T ip T /(r 2 + A 2 ), which approaches the 
bare linewidth of the resonator, k/2, for strong detuning. 

When the coupling is much stronger than the deco- 
herence rate, the linewidth depends nomonotonously on 
the detuning exhibiting an enhancement at the lasing 
transition (green solid line) |11| . Different from the con- 
ventional lasing situation where the emission spectrum 
is mainly broadened by phase fluctuations, in the strong 
coupling regime at the lasing transition, the contribution 
from amplitude fluctuations is magnified via its coupling 
to the phase fluctuations, leading to the peak shown in 
the figure. In this transition regime, the phase correlator 
(a t (t)a(0)) exhibits a non-exponential decay in time. 

A further surprising phenomenon occurs for a strong 
pure dephasing rate r*/2-7r = 4 MHz. Instead of broad- 
ening the emission spectrum, the pure dephasing actu- 
ally reduces the linewidth (red dashed line), leading to 
a much sharper spectrum peak. An intuitive interpreta- 
tion can be obtained by examining the linewidth at reso- 
nance, where a semi-quantum approach is accessible (Tol j . 
In this case, the linewidth at resonance in the deep las- 
ing regime can be estimated as «l oc g 2 /(T v (n)). With 
growing pure dephasing, the photon number changes lit- 
tle while the total decoherence rate T v increases linearly 
with the pure dephasing rate, which leads to the decrease 
in the linewidth. 



V. CORRELATIONS BETWEEN LASING AND 
TRANSPORT PROPERTIES 

A current flows in the system if the tunneling cycle 
|0,0) ->■ |1,0) -> |0,1) -+ 1 0,0) is completed. This can 



be achieved by either coherent interdot tunneling or a 
coherent process involving the excition of photons in the 
resonator. We calculate the current using 



(10) 



where the index i refers to the states \g), \e), and |0,0), 
and denotes the transition rate from state \i) to 

\j). In Fig. [7] we plot the current as a function of the 
dot detuning e. As expected, two resonance peaks show 
up in the current. The coherent interdot tunneling leads 
to the broad peak around e = with width given by the 
tunneling rate t (here t ^> hT) [30|,[3lj]. Interestingly, the 
second peak due to transitions involving the excitation of 
photons is - for realistic values of the parameters - much 
narrower. As shown in Fig. [FJ it correlates with the lasing 
state, since the excitation of a photon in the resonator is 
caused by the electron tunneling between the two dots. 
Both the lasing state and the current peak exist only in 
a narrow "resonance window" |A| < W/2. This window 
can be estimated from Eq. ([7]) with condition (n) a > 0, 
namely, 



W 



32 cos 9 g 2 



[cos(20) + 7]/cr 



-1, 



(11) 



which reduces to W ss 2gy/T/n — 1 for small 9. To esti- 
mate the height of the narrow current peak, we adopt an 
adiabatic approximation assuming the dynamics of the 
resonator to be much slower than that of the quantum 
dots [HI]. For k r and small 9, the peak value of the 
current is given by 



I(A = 0) ~ eTj2 p ( n ) 



2(n+ 1) 



3(n+ 1) +r 2 /(V) 



(12) 



Here P(n) ~ (T / k) P (0)Uf =1 [3l + r 2 /(4 5 2 )]- 1 denotes 
the probability of having n photons in the resonator. 
When the coupling to the resonator is strong compared to 
the incoherent tunneling T, the peak current approaches 
2er/3. 

The correlation between the lasing state and the trans- 
port current is remarkable in two ways. On one hand, the 
current peak, which may be easier to measure than the 
photon state of the resonator, can be used as a probe of 
the lasing state. On the other hand, the rather sharp 
resonance condition needed for the lasing makes the cur- 
rent peak narrow, while at the same time the value of 
the current is reasonably high. This allows resolving in 
an experiment small details of the dot properties. 

The effects of decoherence on the transport current 
are illustrated in Fig. [8j Relaxation, on one hand, re- 
duces the efficiency of the coherent transition between the 
excited and ground states. Hence within the resonance 
window the current induced by the coherent transition 
decreases. On the other hand, relaxation also opens up 
an incoherent channel which increases the current. Pure 
dephasing leads to a broadening of the resonance peak in 
the current, as well as of the photon number. 
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FIG. 8: (Color online) Current as a function of detuning A 
with the green solid line for the result with vanishing deco- 
herence = T* = 0, the blue dotted line with relaxation 
rate Fj./27r = 1.2 MHz but vanishing pure dephasing, and the 
red dashed line with pure dephasing rate T* v /2ti — 4 MHz but 
Ti — 0. The photon number (gray dot-dashed line) is also 
shown for comparison. 

VI. SUMMARY 

We have investigated a double quantum dot coherently 
coupled to a transmission line resonator. A population 
inversion between the dot levels is created by a pump- 
ing process due to a voltage applied across the double 



dot. A lasing state of the radiation field develops within 
a sharp resonance window. It correlates with a peak in 
the transport current. The sharp resonance condition al- 
lows for resolving small differences in the dot properties. 
This opens perspectives for applications of the setup and 
operation principle for high resolution measurements. 

We analyzed the effects of dissipation on the dot lev- 
els. Both relaxation and pure dephasing deteriorate the 
lasing state, accompanied by a decrease in the photon 
number and an increase of amplitude fluctuations. Re- 
laxation processes shorten the life time of the two-level 
system and reduce the population inversion and hence 
the lasing effect. Pure dephasing, which also shortens 
the life time, broadens the effective resonance window be- 
tween the quantum dots and resonator. Surprisingly, the 
emission spectrum can become even sharper with mod- 
erate pure dephasing rate. For the transport current, the 
relaxation opens up an extra incoherent channel which 
increases the current outside the resonance window with 
the resonator. 
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